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We present results from a systematic direct-numerical simulation study of pressure fluctuations 
in an unforced, incompressible, homogeneous, and isotropic, three-dimensional turbulent fluid. At 
cascade completion, isosurfaces of low pressure are found to be organised as slender filaments, 
whereas the predominant isostructures appear sheet-like. We exhibit several new results, including 
plots of probability distributions of the spatial pressure-difference, the pressure-gradient norm, and 
the eigenvalues of the pressure-hessian tensor. Plots of the temporal evolution of the mean pressure- 
gradient norm, and the mean eigenvalues of the pressure-hessian tensor are also exhibited. We 
find the statistically preferred orientations between the eigenvectors of the pressure-hessian tensor, 
the pressure-gradient, the eigenvectors of the strain-rate tensor, the vorticity, and the velocity. 
Statistical properties of the non-local part of the pressure-hessian tensor are also exhibited, for 
the first time. We present numerical tests (in the viscous case) of some conjectures of Ohkitani 
[Phys. Fluids A 5, 2570 (f993)] and Ohkitani and Kishiba [Phys. Fluids 7, 411 (1995)] concerning 
the pressure-hessian and the strain-rate tensors, for the unforced, incompressible, three-dimensional 
Euler equations. 
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I. INTRODUCTION 

The pressure at a point in an incompressible fluid 
is a non-local functional of the velocity field, and in- 
herently difficult to measure in laboratory experiments. 
However, high-resolution direct-numerical simulations 
(DNS) provide accurate statistics of the pressure field. 
Pressure fiuctuations have been extensively studied in 
both numerical studieslll d S H H and laboratory 
experiments 0, S 0, llrf] of statistically steady, homoge- 
neous, and isotropic turbulent flows. High-resolution nu- 
merical studies0^|3, 0, 01 show, that the pressure spec- 
trum exhibits a wavenumber range with power-law scal- 
ing, in accordance with predictions from a phenomeno- 
logical theory due to Kolmogorov|llj. The pressure 
probability distribution is widely accepted to be nega- 
tively skewed, and to exhibit an exponential low-pressure 
tail(3,l3- Regions of low pressure are found to be organ- 
ised as slender filamentary structures in both numerical 
studies 2] and laboratory experiments^. Il0l| . 

The 'canonical' isotropic turbulent flow is a decaying 
turbulent flow behind a grid^J. The study of decay- 
ing turbulence is important since the results are unin- 
fluenced by statistics of the forcing and directly reflect 
effects of the nonlinear terms in the Navier-Stokes equa- 
tions (see Eqs. below). In contrast to statistically 
steady turbulence, systematic numerical studies of pres- 
sure fluctuations within the context of decaying, homo- 
geneous, and isotropic turbulence are extremely scarce. 
The only available work is a low-resolution DNS study 
due to Schumann and Patterson^ who exhibited plots 
of the root-mean-square pressure fluctuations as a func- 



tion of the time, and the isosurfaces of low pressure. The 
low-pressure isosurfaces in this study0] were shown to 
be organised as 'cloud'-like structures, in contrast to the 
slender fllaments seen in DNS studies '5| of statistically 
steady turbulence at higher resolutions. The pressure- 
hessian tensor and the pressure-gradient were not studied 
in this workfT^. In both statistically steady and decay- 
ing turbulence, a comprehensive study of possible align- 
ments between vectors of interest in a turbulent flow, 
namely, the eigenvectors of the pressure-hessian tensor, 
the pressure-gradient, the eigenvectors of the strain-rate 
tensor, the vorticity, and the velocity, is entirely lacking. 
In this paper, we present results from a systematic nu- 
merical study of the pressure, the pressure-gradient, and 
the pressure-hessian tensor in an unforced, incompress- 
ible, homogeneous, and isotropic turbulent fluid. We ex- 
hibit several new results, including plots of the prob- 
ability distributions of the spatial pressure-difference, 
the pressure-gradient norm, and the eigenvalues of the 
pressure-hessian tensor, temporal evolution of the mean 
pressure-hessian eigenvalues and of the mean pressure- 
gradient norm, as well as isosurfaces of the pressure and 
the pressure-gradient norm, at cascade completion. Sta- 
tistical properties of the non-local part of the pressure- 
hessian tensor are also exhibited, for the flrst time. 
We construct the general alignment picture between the 
eigenvectors of the pressure-hessian tensor, the pressure- 
gradient, the eigenvectors of the strain-rate tensor, the 
vorticity, and the velocity. Ohkitanij^^ ^ind Ohkitani 
and Kishiba 13 have derived several interesting results 
for the unforced, incompressible, three-dimensional, in- 
viscid Navier-Stokes equations (the Euler equations) con- 
cerning the pressure-hessian and the strain-rate tensors. 
We exhibit numerical tests of some conjectures for the 
Navier-Stokes case. 
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The unforced Navier-Stokes equations are 



dv 

'dt 



+ (vV)v = -^ + z.V2v, (1) 
P 



where v is the kinematic viscosity, p is the (constant) 
density, and p is the pressure. On taking the divergence 
of Eqs. (0, and enforcing incompressibihty V ■ v — 0, 
the pressure is determined by 
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(2) 



where the enstrophy = UiUJi, tUi = CijkdjVk is the 
vorticity (e^fe is the Levi-Civita tensor), i,j,k = 1,2,3, 
with a summation imphcit over repeated indices, and 
= SijSij, Sij = \/2{djVi + diVj) is the strain-rate 
tensor. 



II. NUMERICAL METHOD 

We use a pseudospectral method to solve Eqs. ^ 
numerically, in a cubical box of side 27r with periodic 
boundary conditions and 224'^ Fourier modes. In this 
paper, we do not address issues pertaining to the scal- 
ing of higher-order structure functions of the pressure- 
difference (or pressure-velocity correlations) or investi- 
gate dissipation-scale properties, and believe that our 
spectral resolution is adequate for the types of studies 
that we have carried out (barring the pressure spec- 
trum, see below). We have checked that our results 
are unaffected by resolution considerations, by compar- 
ing with results from a 128'^ DNS study with identi- 
cal (initial) Reynolds number. For the temporal evo- 
lution, we use an Adams-Bashforth scheme (step size 
5t = 10~^) with double-precision arithmetic and set 
p = 1, V — 10"^. We include a hyperviscous term 
of the form VhV'^'v in Eqs. iQ), with Vh = 10~^ and 
have explicitly checked that our results are unaffected 
by the inclusion of hyperviscosity. We note that Borue 
and Orszag0 have carried out a 256'^ DNS study of 
decaying, isotropic turbulence with hyperviscosity, and 
they conclude that "...inertial-range dynamics may be 
independent of the particular mechanism of small-scale 
dissipation..." (p. R859, Ref. ^3). The initial velocity 
field is taken to be v(k, <o) ^ k'^e~^ e'-^^ {k = |k| is the 
wavenumber), with 0k independent random variables dis- 
tributed uniformly between and 27r. This corresponds 
to an initial kinetic-energy spectrum E{k,to) ~ k'^e~'^^ 
(with E{k, t) = |v(k, i)p the one-dimensional spectrum), 
which is a convenient choice that develops a cascade to 
large wavenumbers (see below). We measure time in 
units of the initial 'box-size' time tq = 27r/wJ^„j, (here 
To equals 4.02), v°^, = [(X;^ \'^{^M)?)Y''^ is the root- 
mean-square value of the initial velocity, with the dimen- 
sionless time t = t/To (t is the product of the number of 
steps and 5t). We define Re^ = 2t:v^^^Iv to be the value 
of the initial 'box-size' Reynolds number (here i?eo equals 



982464). Our results are obtained for times to < t << t*, 
where t^, is the time at which the (growing) integral 
scale L{i) = ((Ek |v(k, t)P A)/ Ek |v(k, i)P) becomes 
of the order of the linear size of the simulation box. For 
times t >t^,, finite-size effects which might well be non- 
universal, modify the numerical results, and are not con- 
sidered here. 

In Fig. ^ we show some preliminary results that serve 
as a check of our numerical method and parameter val- 
ues (which were chosen to ensure linear stability of the 
numerical scheme). Figure Qfa) shows on a log-log plot, 
the scaled kinetic energy spectrum k^/^E{k, r) as a func- 
tion of the wavenumber k. On starting with the spec- 
trum specified above, a cascade of energy is seen to large 
wavenumbers. The plots are equispaced in time with a 
temporal separation of t = 0.24. The plot with open 
circles is calculated at cascade completion at the dimen- 
sionless time r = Tc = 0.71, and shows a wavenumber 
range (for 1 < k 5. 10) that exhibits the well-known 
—5/3 power-law [nl [13, [13. Upon cascade completion, 
the shape of the energy spectrum does not change ap- 
preciably (except at large wavenumbers where it falls), 
but the kinetic energy decays monotonically. In Fig. 
n^b), we plot the normalized kinetic energy-dissipation 
rate e(r)/eo [e(t) = Ek ^^l'^(kj *)P] as a function of the 
dimensionless time t. The kinetic energy-dissipation rate 
peaks 13 01 at r = Tc, corresponding to cascade comple- 
tion in the kinetic energy spectrum, and decreases there- 
after. The turbulence may be considered as 'fully devel- 
oped' at T = Tc and our spatial results (see below) will 
be calculated at this instant of time. 



III. NUMERICAL RESULTS 

A. Pressure 

At each grid point, we compute p by Fourier- 
transforming Eq. |(2J), solving for the pressure, and 
inverse Fourier-transforming to physical space. Kol- 
mogorov phenomenology predicts that the pressure 
spectrum in statistically steady turbulence exhibits a 
wavenumber range with the power-law scaling p{k) ~ 
/c~^/''|^. Ishihara, et.aZ.Q have confirmed the scal- 
ing law in a 2048^^ DNS study of statistically steady, 
homogeneous, and isotropic turbulence at a Taylor- 
scale Reynolds number Re\ — 732, whereas Tsuji and 
Ishihara[l3 have observed the scaling law by measuring 
pressure fluctuations in the centre line of a freely decaying 
turbulent jet. In Fig. |2Ia), we plot the scaled pressure 
spectrum k'^/^p{k, t) as a function of the wavenumber fc, 
at cascade completion. In order to observe Kolmogorov- 
type scaling in the pressure spectrum over a substan- 
tial wavenumber range, a considerably higher (initial) 
Reynolds number is reauired|20j as compared to the re- 
sult for the kinetic-energy spectrum [cf. plot with open 
circles in Fig. ^a)]. The spectral resolution of our study 
is inadequate for the purposes of fitting a power-law, and 
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FIG. 1: (a) Log-log plot of the temporal evolution of the 
scaled kinetic energy spectrum k^^^E{k, r) as a function of 
the wavenumber k at temporal separations of r = 0.24. The 
plot with open circles is calculated at cascade completion, at 
dimensionless time r = Tc = 0.71. 

(b) Plot of the normalized kinetic energy-dissipation rate 
e(r)/eo as a function of the dimensionless time r. 



the pressure spectrum is found to exhibit a —7/3 power- 
law only in the narrow wavenumber range 20 < fc < 50. 

In Fig. |2Ib), we plot the normalized probability distri- 
bution 'P(Apr) of the pressure-difference Ap^ = p(x -I- 
r) — p(x) at cascade completion, for grid-spacing val- 
ues r = |r| = 1,50. For the large seperation r = 50, 
P(Ap5o) is found to be close to a Gaussian distribution 
[cf. dashed-line curve in Fig. |2Ib)] with a skewness equal 
to zero, and a kurtosis equal to 3.49. For r = 1, stretched- 
exponential tails which are roughly symmetrically placed 
about Api — are observed, with 7^(Api) having a kur- 
tosis equal to 7.80. We note that 7'(Api) does not ex- 
hibit a Gaussian core at small values of Api/cr {a is the 
standard deviation). Our results for V{Apr) resemble 
those obtained for probability distribution of the spatial 
velocity-difference|2l| at large and small grid-spacings 
respectively. Cao, e<.aZ.|3| obtained similar results for 
■p(Apr) from a 512'^ DNS study of statistically steady, 
homogeneous, and isotropic turbulence at Re\ = 218. 

In Fig. Eta), we plot the normalized probability distri- 
bution Vip) of p, at cascade completion. The mean pres- 



sure (p) (angular brackets denote a volume average) was 
found to equal zero, as is expected [Tlj| for isotropic tur- 
bulence, whereas the skewness was found to equal —0.11 
and the kurtosis equalled 5.62. Holzer and Siggia'^^] have 
shown analytically, that the probability distribution of 
the pressure is negatively skewed and has an exponential 
tail, for a Gaussia'n}2^ velocity probability distribution. 
Brachet,2^ has found a low-pressure exponential tail in 
an 864^ DNS study of decaying, isotropic turbulence 
with Taylor- Green |2^ initial conditions. The error bars 
in this studyli^l are probably larger than those in the 
data of Fig. Ola). Both PumirQ and Cao, ei.aL0] ob- 
served a stretched-exponential tail at low pressures, and 
a roughly Gaussian tail at high pressures, in DNS stud- 
ies of p in statistically steady turbulence. We confirm the 
result for negative pressures in the case of decaying tur- 
bulence, where the fit T'dp]) - e^'^lpl", /3 = 2.55 ±0.01, 
a = 1.35 ± 0.01 (error-bars from least-square fits), is ob- 
served at cascade completion. However, at positive pres- 
sures, a stretched-exponential tail with /3 = 1.80 ± 0.01, 
a — 1.52 ± 0.02 (error-bars from least-square fits) is ob- 
served in our study. 

In Fig. |3Ib), we plot iso-p surfaces for the isovalue 
p = {p) at cascade completion, which appear to be crum- 
pled sheet-like structures (found throughout the isovalue 
range [{p) — cr, (p) + cr]). Equation ^ suggests that re- 
gions with high vorticity and low strain-rates are simul- 
taneously sources of low pressure. Such regions with in- 
tense vorticity have been observed by Douady, ei.aZ.|3| 
and Villermaux, etaZ. jTfll| in statistically steady turbu- 
lence experiments, by using cavitation as a visualization 
technique, in a liquid seeded with bubbles. Schumann 
and Patterson|l2j exhibited plots of low-pressure isosur- 
faces from a 32^^ DNS study of the unforced, incompress- 
ible, three-dimensional Navier-Stokes equations, which 
were shown to be organised as 'cloud'-like structures. In 
our study, at early times r << Tc, regions of low pres- 
sure (with the isovalue p ^ (p) ^ 2(t) are found to be 
sheet-like [see Fig. EJa)]. At cascade completion, low- 
pressure regions are found to be organised as slender fil- 
aments [see Fig. E^b)], with diameter of the order of the 
grid spacing, and a contour length that occasionally ex- 
tends nearly to the linear size of the simulation box. We 
choose to quote dimensions of the structures relative to 
the (fixed) box-size and the grid-spacing, since both the 
Kolmogorov (dissipative) and the integral length scales 
vary in time, in decaying turbulence. Isosurface plots 
of the enstrophy and the squared strain-rates show (see 
Ref. 0) that highly-strained regions occur close to re- 
gions of intense enstrophy, and Eq. ||2Jl suggests a lack of 
any well-defined fluid-mechanical structure of the high- 
pressure regions. In our study, iso-p surfaces in the range 
p > {(p) + a) were indeed found not to exhibit any par- 
ticular structure at cascade completion. 
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FIG. 2: (a) Log-log plot of the scaled pressure spectrum 
k^^^pik, t) as a function of the wavenumber k, at cascade 
completion. The horizontal line is drawn for reference, 
(b) Semilog plot of the normalized probability distribu- 
tion ViApr) (cr is the standard deviation) of the pressure- 
difference Apr at cascade completion, for grid-spacing values 
r = 1, 50. The dashed- line curve is a normalized Gaussian 
distribution for comparison. 

B. Pressure-Hessian Tensor 

The pressure-hessian tensor Pij = dijp appears [TslIT^ 
in the evolution equation for the strain-rate tensor Stj. 
At each grid point, we compute the eigenvalues Ai.p, 
A2.P, and Aa^p (with the convention Ai_p > A2.P > Aa^p) 
of Pij and the corresponding orthonormal eigenvectors 
/i, /2, and /a. We also compute the eigenvalues Ai,^, 
A2.S', and As^s (with ordering from extensional to com- 
pressive strain-rates Ai^s > A2,s > ^3,s) of Sij and the 
corresponding orthonormal eigenvectors ei, 62, and 63. 

In Fig. Efa), we plot the normalized probability distri- 
bution V{Xi^p) of the eigenvalues A^^p, at cascade com- 
pletion. The skewnesses of the eigenvalue distributions 
were found to equal 4.73, 5.25, and —4.13 for i = 1,2, 
and 3 respectively. 

In a constant-density flow, incompressibility requires 
that '^iXi,s = 0; however, there is no such constraint on 
Ai,p. The inset plot in Fig. Efb) is the normalized prob- 
ability distribution of the trace of Pij, V[x — Tri^Pij)] 
{Tr{Pij) = J2k-^k,p) at cascade completion, which is 
roughly symmetrically placed about x ~ 0. As is well- 
known in both decavingfisj and statistically steadv[2H 
turbulence, regions of intense vorticity (say, for iso-|u;| 




FIG. 3: (a) Semilog plot of the normalized probability dis- 
tribution 'P(p) of the pressure p, at cascade completion. The 
dashed-line curve is a normalized Gaussian distribution for 
comparison. 

(b) Plot of iso-p surfaces for the isovalue p = (p), at cascade 
completion. 



values greater than (|a;|) + 2a) are found to be organ- 
ised as filamentary structures, and are spatially more lo- 
calised than regions of high strain-rate. Equation Q 
suggests that locally, in regions of intense enstrophy, 
Tr{Pij) = V^p > 0. In Fig. Etb), we plot the normalized 
probability distribution V\Tr{Pij)] of the trace of Pij at 
cascade completion conditioned on > (|w|)-t-2(T, which 
is found to have a positive mean as expected. 

Ohkitani and Kishiba^^, have shown that the 
pressure-hessian tensor Pij can be decomposed into the 
sum of a diagonal tensor 5ij\/^p/3 (the 'local' term), 6ij 
is the Kronecker delta, and a symmetric, zero-diagonal 
tensor (5ij^3 (the 'non-local' term). The local term, 
as the name indicates, can be expressed purely in terms 
of the vorticity and the strain-rate at each point of the 
fluid [see Eq. 0], however, the non-local term can be 
expressed only in terms of an integral over the entire 
fluid volume. Ohkitani and KishibaJ4j have also shown 
(for Tavlor-Green[25j initial conditions) that the non- 
local term contributes significantly to enstrophy growth. 
Since Qij — Pij — iJ^ V^p/3 has zero trace, its eigenvalue 
Ai,Q > 0, Aa^Q < 0, and the sign of A2,q is indeter- 
minate. In Fig. Efc), we plot the normalized probabil- 
ity distribution 'P{\i.Q) of the eigenvalues Ai^g, at cas- 
cade completion. We find that A2,q has a positive mean. 
The statistically preferred ratio of the mean eigenvalues 
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FIG. 4: (a) Plot of iso-p surfaces for the isovalue p = (p) — 2a 
at the dimensionless time r << Tc- 

(b) Plot of iso-p surfaces for the isovalue p — (p) — 2a, at 
cascade completion. 



{Xi,q) ■ {^2,q) '■ {^3,q) was found to equal 46 : 1 : —47, 
at cascade completion p^. 

In Fig. inia), we plot the mean eigenvalues {Xi,p) as 
a function of the dimensionless time r, these evolve in 
a way that is similar to the temporal evolution of the 
kinetic-energy dissipation rate, with a peak in the mag- 
nitude of (Ai^p), at cascade completion [cf. Fig. H^b)]. 
Similar results (not shown here) were obtained for the 
temporal evolution of (Ai^g). Ohkitani[13i|. in a 128'^ DNS 
study of the unforced, incompressible, three-dimensional 
Euler equations, showed that A3.P changes sign (locally, 
in regions of intense enstrophy) from positive to nega- 
tive at early times. However, in our study, we find that 
(A2,p) changes sign [see Fig. Elb)] at r = 0.37 from neg- 
ative to positive, whereas (Aa^p) remains negative at all 
times (not shown here). 

In Fig. CJa), we plot the normalized probability dis- 
tribution of the cosines of the angles between the eigen- 
vectors fi of Pij^2^ and the pressure-gradient Vp (see 
below), at cascade completion. In Fig. [7Jb), we plot the 
normalized probability distribution of the cosines of the 
angles between the eigenvectors fi and uj, at cascade com- 
pletion. Ohkitani 13.] showed that lo is preferentially par- 
allel (or anti-parallel) to the eigenvector /a corresponding 
to the pressure-hessian eigenvalue As^p smallest in mag- 
nitude, in contrast to our result, which shows that uj is 
preferentially parallel (or anti-parallel) with the eigenvec- 
tor /2 corresponding to the intermediate eigenvalue A2,p. 
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FIG. 5: (a) Semilog plot of the normalized probability distri- 
bution V{Xi,p) of the eigenvalues Xi,p of the pressure-hessian 
tensor Pij, at cascade completion. 

(b) Semilog plot of the normalized probability distribution 
V[x = Tr(P„)] of the trace of Pij {Tr{P,j) = Efe Afc.p) at 
cascade completion, conditioned on |a;| > -|- 2a. The in- 
set is a semilog plot of the normalized probability distribution 
P{x) of a; = Tr{Pij), at cascade completion. 

(c) Semilog plot of the normalized probability distribution 
P{Xi,Q) of the eigenvalues Ai,Q of the tensor Qij (see the text 
for the definition), at cascade completion. 



In Fig. [7fc), we plot the normalized probability distribu- 
tion of the cosines of the angles between the eigenvectors 
fi and the velocity v, at cascade completion. 



C. Pressure Gradient 

A systematic numerical study of the pressure gradient 
is entirely lacking in both statistically steady and decay- 
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FIG. 6: (a) Plot of the mean eigenvalues (K^p) of Pij as a 
function of the dimensionless time r. 

(b) Plot of (A2,p) as a function of the dimensionless time r, 
in the range < r < 0.75. 



we plot the normalized 
of the Euclidean norm 
for vector x with components Xi) of the 



ing turbulence. In Fig. |HIa), 
probability distribution P(|Vp 

(|x| ^ VE:^ 
pressure-gradient Vp at cascade completion, which does 
not exhibit a stretched-exponential tail. The author has 
been unable to determine a functional form which gives 
a good fit for the tail of T'dVp]). 

In Fig. IHIb), we plot the mean pressure-gradient norm 
(|Vp|) as a function of the dimensionless time r, which 
is observed to peak at r = Tc, as does the kinetic-energy 
dissipation rate [see Fig. Qlb)]. 

In Fig. IHIc), we plot iso-|Vp| surfaces for the isovalue 
\Vp\ = (|Vp|) -I- 2a, at cascade completion. The iso- 
surfaces of intense pressure-gradient, which are filamen- 
tary in shape, are found to resemble the /o^/^pressure 
isosurfaces in Fig. Elb). Iso-|Vp| surfaces in the range 
Vp| < ((|Vp|) — a) were not found to exhibit any par- 
ticular structure at cascade completion. 

In Fig. Eta), we plot the normalized probability dis- 
tribution of the cosines of the angles between Vp and 
the eigenvectors of Sij. Ashurst, et.al. S], noted a ten- 
dency for the alignment of Vp (caused by velocity fluctu- 
ations alone) with the most compressive strain direction 
63 in a 128'^ DNS study of a statistically steady turbu- 
lent shear-flow. However, in our study, Vp is not found to 
be preferentially parallel (or anti-parallel) with 63, and 
we observe a peak in the magnitude of 7'[cos(Vp, 63)] 
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FIG. 7: (a) Plot of the normalized probability distribution of 
the cosines of the angles between the eigenvectors fi of Pij 
and the pressure gradient Vp, at cascade completion. 

(b) Plot of the normalized probability distribution of the 
cosines of the angles between the eigenvectors fi and the vor- 
ticity uj, at cascade completion. 

(c) Plot of the normalized probability distribution of the 
cosines of the angles between the eigenvectors fi and the ve- 
locity v, at cascade completion. 



at |cos(Vp, 63)1 « 0.86, which indicates a preferential 
relative angle « 7r/6. In Fig. El^b), we plot the normal- 
ized probability distribution of the cosine of the angle be- 
tween Vp and w, at cascade completion. In Fig. EJc), we 
plot the normalized probability distribution of the cosine 
of the angle between Vp and v, at cascade completion. 
Both UJ and v are observed to be preferentially perpen- 
dicular to Vp, in agreement with corresponding results 
from a 128'^ DNS study of the unforced, incompressible, 
three-dimensional Euler equations, due to Ohkitani and 
KishibaQ. 
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FIG. 8: (a) Semilog plot of the normalized probability distri- 
bution ■p(IVpl) (|Vp| is the Euclidean norm of the pressure- 
gradient), at cascade completion. The dashed-line curve is a 
normalized Gaussian distribution for comparison. 

(b) Plot of the mean pressure-gradient norm (|Vp|) as a func- 
tion of the dimensionless time r. 

(c) Plot of iso-|Vp| surfaces for the isovalue |Vp| = {\Vp\)+2a, 
at cascade completion. 
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FIG. 9: (a) Plot of the normalized probability distribution of 
the cosines of the angles between Vp and the eigenvectors 
of the strain-rate tensor Sij, at cascade completion. 

(b) Plot of the normalized probability distribution of the co- 
sine of the angle between Vp and uj, at cascade completion. 

(c) Plot of the normalized probability distribution of the co- 
sine of the angle between Vp and v, at cascade completion. 



D. General Alignment Picture 

Ohkitani and Kishiba^J have found the orientations 
(locally, in regions of intense enstrophy) amongst the set 
of vectors (e^, excluding the pressure gradient and 

the velocity. In Figs. [7|and|51 we have plotted possible 
alignments between the sets of vectors (/i, Vp) and (e^, 
V, Lu) and for completeness, we plot in Figs. 1101 and [TTl 
remaining alignments amongst the set of vectors (e^, v, 
uj) and amongst the eigenvector bases (fi, e^), at cascade 
completion. 

As is well-known, in both decaying^^ and statisti- 



cally steady'2ll turbulence, there is an increased proba- 
bility for alignment (or ant i- alignment) of the interme- 
diate strain-rate eigenvector 62 with lo, relative to the 
alignments between ei and 63 with uj. In Fig. Iior a). we 
plot the normalized probability distribution of the cosines 
of the angles between to and at cascade completion, 
which reaffirms this result for decaying turbulence. In 
Fig. llOf b). we plot the normalized probability distribu- 
tion of the cosines of the angles between v and e^, at 
cascade completion. In Fig. llOf c). we plot the normal- 
ized probability distribution of the cosine of the angle be- 
tween V and w, at cascade completion. It is interesting to 
note that the alignment plots in Figs. [7||51 and liur a).fb) 
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are roughly symmetrically placed about cos6' = 0, how- 
ever Fig. Iiur c) is distinctly asymmetric, with a greater 
probability for v and u to be anii-parallel. This asym- 
metry has been noted in a laboratory experiment [29j of 
decaying turbulent flow past a grid, and is plausibly due 
to effects of the kinetic helicity 30] (an invariant of the 
Euler equations) on the decay process. 

In Fig. [TlT a). 7'[cos(/i, ei)] is found to peak at 
|cos(/i,ei)| 0.71, which indicates a preferential rel- 
ative angle ~ 7r/4, in agreement with corresponding re- 
sults due to Ohkitani and KishibajQ]. The only distinct 
feature in Figs. Illf b'l.fc'l is that /2 is preferentially par- 
allel (or anti-parallel) to 62 (and perpendicular to 63). 

Ohkitanifi^ has conjectured that the pressure-hessian 
tensor P (with components Pij) and the strain-rate ten- 
sor S (with components Sij) in general are not com- 
mutative, and therefore cannot be simultaneously diag- 
onalized. It is of interest to determine the (statistically 
preferred) relative orientation, between the two frames 
with respect to which S and P are diagonalized. Ohki- 
tani and Kishiba[T^ have conjectured that the configura- 
tion of relative alignments between the eigenvector bases 
{fii^i) is one with "least commutativity between S and 
P out of all possibilities with one axis in common" (p. 
414, Ref. [IJI) at cascade completion. In order to test 
these conjectures, we choose the standard matrix norm 
\\A\\ = (maximum eigenvalue of A"^ A)^^"^ ^1^, where the 
superscript T denotes the transpose-conjugate. In Fig. 
1121 we plot the mean norm (nils', P]||) of the commutator 
[S, P] — SP — PS* 32] as a function of the dimension- 
less time T and find that the value peaks at t = 0.71. 
The value of the mean norm depends on the choice of 
the norm; however, the trend (which is independent of 
this choice) indicates that the relative configuration of 
the eigenvector bases (/i,ei) is "least commutative" at 
T — 0.71, which is equal to the time at which the ki- 
netic energy-dissipation rate peaks, in accord with the 
Ohkitani and Kishiba^^ conjecture. 



IV. CONCLUSION 

To summarize, we have presented results from a sys- 
tematic numerical study of pressure fluctuations in an un- 
forced, incompressible, homogeneous, and isotropic tur- 
bulent fluid. At cascade completion, isosurfaces of low 
pressure are found to be organised as slender filaments, 
whereas the predominant pressure isostructures appear 
sheet-like. We have exhibited several new results, in- 
cluding plots of the probability distributions of the spa- 
tial pressure-difference, the pressure-gradient norm, and 
the eigenvalues of the pressure-hessian tensor, at cas- 
cade completion. Plots of the temporal evolution of the 
mean pressure-gradient norm, and the mean eigenvalues 
of the pressure-hessian tensor have also been exhibited. 
We have found the statistically preferred orientations be- 
tween the eigenvectors of the pressure-hessian tensor, the 
pressure- gradient, the eigenvectors of the strain-rate ten- 




cos ( V, ft) ) 

FIG. 10: (a) Plot of the normalized probability distribution 
of the cosines of the angles between uj and the eigenvectors 
ei, at cascade completion. 

(b) Plot of the normalized probability distribution of the 
cosines of the angles between v and the eigenvectors ei, at 
cascade completion. 

(c) Plot of the normalized probability distribution of the co- 
sine of the angle between v and w, at cascade completion. 



sor, the vorticity, and the velocity at cascade comple- 
tion. Statistical properties of the non-local part of the 
pressure-hessian tensor have also been exhibited, for the 
first time. We have presented numerical tests (in the vis- 
cous case) of some conjectures for the unforced, incom- 
pressible, three-dimensional Euler equations, proposed in 
earlier studies. 
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FIG. 11: (a) Plot of the normalized probability distribution 

of the cosines of the angles between the eigenvectors /, and 
ei, at cascade completion. 

(b) Plot of the normalized probability distribution of the 
cosines of the angles between the eigenvectors fi and 62, at 
cascade completion. 

(c) Plot of the normalized probability distribution of the 
cosines of the angles between the eigenvectors fi and 63, at 
cascade completion. 




FIG. 12: Plot of the mean norm (||[S',P]||) (see the text for the 
definition of the norm) of the commutator [S, P] = SP — PS 

of the strain-rate S and the pressure-hessian P tensors, as a 
function of the dimensionless time t. 
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